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Abstract 

We study the fermion sign problem in a theory of non-relativistic fermions with a spin- 
independent repulsive interaction. We work in polar co-ordinates in momentum space, which 
makes it straightforward to keep only the low-energy degrees of freedom close to the Fermi surface. 
This is sufficient for the purpose of calculating many physically important low-energy observables. 
We find indications that the sign problem in this effective theory will be weaker than in the full 
theory, so low-energy properties of the theory could be calculated by modifying the action to make 
it positive semi-definite and including reweighting factors in the observables. We discuss suitable 
modifications of the action, and describe a possible lattice realization of the polar momentum space 
formulation of the theory. 
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I. INTRODUCTION 



One of the few known first-principles techniques for calculating the properties of systems 
of strongly-interacting fermions is numerical evaluation of the functional integral. The stan- 
dard approach, which involves Monte Carlo importance sampling of field configurations on 
a space-time lattice, has been very successful in zero-density and high-temperature systems, 
but at non-zero fermion density and low temperature it becomes impractical because the 
fermionic part of the integrand is not positive and therefore cannot be treated as a weight 
for the importance sampling (see, for example, [H, S 0] ) • This has become known as the 
fermion sign problem. 

In this article we investigate the feasibility of numerically calculating low-energy proper- 
ties of a degenerate Fermi system by focusing on fermionic degrees of freedom close to the 
(perturbative) Fermi surface [sl, 0] . This approach, known as Landau Fermi liquid effective 
theory, allows one to study low-energy properties of the degenerate Fermi system, such as 
heat capacity, density response, and transport coefficients, which are dominated by the low- 
energy degrees of freedom (see, e.g., 0]). We expect this approach to be valid as long as the 
coupling is not so strong as to cause major disruption of the Fermi surface, introducing sen- 
sitivity to degrees of freedom far from the perturbative Fermi surface. As long as this is the 
case, the fermion modes far from the Fermi surface are represented by irrelevant operators 
(see, e.g., jsl). Some observables that are sensitive to the degrees of freedom far from the 
Fermi surface, cannot be studied in this approach; these include bulk properties like fermion 
density and energy density. 

We will work in spherical co-ordinates in momentum space. This makes it straightforward 
for us to write down the leading order terms of the Landau Fermi liquid effective theory, 
which includes only those degrees of freedom that have momenta close to the perturbative 
Fermi surface, and to distinguish between interactions that change radial momentum p 
(distance from the Fermi surface, which is defined by the Fermi momentum pp) and those 
that change the angular component. The price we pay is that the interaction term in the 
effective lagrangian is then non-local. It has been argued by Hong and Hsu [sl, [6| that 
the corresponding effective theory for QCD has a positive semidefinite fermion determinant 
to leading order in 1/pf, if gluon interactions that change a fermion's momentum in the 
angular direction are suppressed. The theory we study here is not QCD: it consists of non- 
relativistic fermions interacting via an auxiliary scalar field. However, it offers a convenient 
toy model for using the Landau Fermi liquid effective theory to study the sign problem in 
dense matter, and it also has physical relevance in its own right. We argue that for this 
theory, without any suppression of scattering between angular patches, the effective theory 
will have a more moderate sign problem than the full theory has. Then relatively minor 
changes to the action (discussed in Sec. IIII|) will yield a positive-definite action that could 
be used to generate an ensemble, and observables in the original theory could be calculated 
by reweighting. We hope that future work will investigate this possibility. 

Momentum lattice approaches have previously been applied to gauge theories jo]. 



111 . |12| . but we will restrict our treatment to non-relativistic fermions that feel a spin- 
independent repulsive interaction, mediated by an auxiliary bosonic field A. In Sec. [TTl we 
write down the fermion matrices in momentum space, in polar co-ordinates, for such a 
theory, and show how the fermion sign problem manifests itself. In Sec. IIIII we discuss ways 
to construct a modified theory with a positive semidefinite fermion matrix. In Sec. IIVI we 
outline a lattice realization of the polar momentum coordinate formulation of the theory. A 
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more detailed discussion of the lattice formulation is given in Appendix |X1 



II. THE NON-RELATIVISTIC FERMION ACTION 
A. Fermions in a background scalar field 

The Euclidean time action for a non-relativistic spin-i fermion ip of mass m in d spatial 
dimensions, at temperature T = 1/(3, interacting with a background configuration A[x) 
(where x = (t, x)) of an auxiliary real bosonic field, is 

'S'fcrm = I dt 







d'^x^pUdr ^l + lAix)]'^|J (1) 

V 2m / 



The field has spin-up and spin-down components, but the action is spin-independent so 
we suppress the spin indices. The bosonic field A has its own action Siaix{A) which depends 
on details of the interaction between the fermions. Its specific form will not concern us 
here. If Eq. ([1]) were obtained by a Hubbard- Stratonovich transformation from a purely 
fermionic action with a 4-fermion repulsive interaction, then the specific form of the 4- 
fermion interaction would determine S ^^nj A). In the simplest case the A field could be a 
discrete spin variable with Sauxl^) = [l3|. 

We now rewrite this action in terms of a Nambu-Gor'kov-like field \E', defined by '^^{x) = 
ip-\{x), "^[{x) = ipl{x). The fermion matrix is still diagonal in spin space but no longer 
proportional to the identity matrix, so we display it explicitly: 







Sfrru, = I dt I d'^X^^ 



+ £ + - iA{x 



* (2) 



The partition function of the theory is then 

Z = JdA detMfe,^(A) exp(-5,,.(A)). (3) 

where Mferm is the fermion matrix, given in square brackets in ([2]). 

Because we assumed the repulsive interaction to be spin-independent, Mferm can be split 
into two independent diagonal blocks, M| and Mi, for spin up and for spin down respectively, 
whose determinants are, for general A{x) and /i (including /i = 0) complex. Therefore the 
partition function ([3]) has a sign problem, because the determinant of the fermion operator 
is the square of a complex number, which is in general not a positive number, and therefore 
cannot be used as a Monte-Carlo importance weight. However, when discretized on a spatial 
lattice of spacing a there is a special value, yUpos = d/{ma^), at which this model has no sign 
problem because each spin component has a real fermion determinant. At /i = /ipos the 
diagonal part of the lattice operator is cancelled, leaving only the hopping terms. Then, 



as in the half-filled Hubbard model [Ij] , there is a particle-hole symmetry, and an associated 
transformation on the fermions that relates M| to its complex conjugate, and similarly for 
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B. Fermions in polar momentum space 



We now Fourier transform in space (not in time), 

^(x) = / e^P-vl/(t,p), A{x) = [ e^P-M(t,p), [ ^ [ ^^^^^ 
Jp J p Jp J 



(27r)' 



(4) 



where the momentum-space field has components ^'(t,p) = {ip^(t,p),ip'^^(t, — p)), obtaining 

Siena = P) ^* + 1^ " ^''(P " + P ~ "l)) ^t(^. q) 



+ ^i(^,P) 



;5 P^ 



5'(p - q) - -P + q)^ ^{(t,q) 



(5) 



Now let us go to spherical coordinates in the spatial momentum. We write each momentum 
in terms of a direction v and a distance p from the perturbative Fermi surface, defined by 
2m/i, 



P = {Pf + p)v, 



+A 



{pF+pY~^dp 



(6) 



-A 



where the factor of l/(27r)'^ is included in the measure of the angular integral. One could 
think of the angular variable v as indexing "patches" on the Fermi surface. We have imposed 
an ultraviolet cutoff A ^ p^? on p, so only modes close to the Fermi surface are included. To 
obtain physical predictions one must integrate out modes with \p\ > A in the standard way 
(see, e.g., Ref. [Sl]). Many of the low-energy observables of interest to us, such as transport 
properties, vanish in the T — >■ limit and are therefore insensitive to the U.V. cutoff (recall 
that there are no temperature-dependent divergences in a quantum field theory) as long as 
A is much greater than their characteristic momentum scale, typically of order the thermal 
momentum mT /pp. 

Our momentum-space approach will therefore be valid as long as fermionic degrees of 
freedom near the Fermi surface play the dominant role in low energy physics. At sufficiently 
strong coupling it may become invalid, for example if there were Bose-Einstein condensation 
of spatially bound two-fermion states. Our approach is therefore most relevant to the region 
of parameter space where the coupling is too strong to allow perturbative approaches, but 
not so strong as to bring in degrees of freedom far from the Fermi surface. An example 
would be the onset of superfiuidity in a moderately-strongly-coupled Nambu-Jona-Lasinio 



model [15|]. It is important to note that the momentum-space approach can be used to map 
its own region of validity. As long as there is a range of values of the cutoff A, obeying 
A <^ pp., where calculated values of low-energy observables such as transport properties are 
independent of the cutoff, we can conclude that fermion modes near the Fermi surface are 
dominating the physics. If we find that we have to push the cutoff up to values of order pp, 
and there is still A-dependence, then our assumption has broken down. 
The fermion action ^ becomes 



S^ 



fcrm 



/ / dpdq {pp + pY'^ipp + qf-^ ^\t, u,p) 

J t,u,v J 



M| 




^{t,v,q) 



(7) 
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where '^(t,v,q) = {tpi{t,v,q),ip^^{t,v,q)), and 



M- 



M 



dt 



Pf_ 


p + 


p2 


m 


2m 


Pf_ 


P- 


p2 


m 


2m 



Su,vSp^q + iA{t, pf{u -v) +pu- qv) 
Su,vSp^q - iA{t, -pf{u -v)-pu + qv) . 



(8) 



where Su^^Sp^g is a simphfied notation for 6'^{p — q) written in spherical coordinates. M| 
describes particles with spin up, and describes holes with spin down. We have explicitly 
shown the indices to make it clear that each is a matrix in momentum space. The indices 
are u and v, which specify directions in momentum space, and p and q, which specify the 
distance from the perturbative Fermi surface. 



C. Fermion determinant for low-energy modes 

We now focus on the degrees of the freedom near the Fermi surface, which we expect to 
be most important for the low-energy properties of the system. To do that, let us represent 
(IH]) and the integration measure ([6]) as an expansion in p/pr, where p is the distance in 
momentum space from the perturbative Fermi surface. This will eventually become an 
expansion in mT/p'jp (equivalently T//i) for observables that are not sensitive to the cutoff, 
or in A/pf where A is the ultraviolet cutoff on p. In the rest of this paper we will focus on 
the leading order contribution. The measure then becomes 

/ =P?"'^ / rdp, (9) 

Jp J V J —A 

and the fermion matrices are 

Su,vSp,q + iA(t^ Pf{u -v) +pu- qv) , 

h,vSp,q - iA{t, -Pf{u -v)-pu + qv) . (10) 

Eqs. (l9i).( fT0l) are obtained by integrating out degrees of freedom with \p\ > A. This means 
that irrelevant couplings (higher-dimension operators, such as higher derivative terms) are 
induced, and relevant couplings such as the fermion mass (and also the parameters of the 
bosonic field action S'aux) must be modified to follow a line of constant physics. We have 
dropped the irrelevant terms since they are suppressed by powers of A/Aq where Aq is the 
cutoff in the original theory, and in (fTOj) we have written the mass as m(A) as a reminder 
that m flows with the cutoff. In the remainder of this paper we will not show the cutoff- 
dependence of the relevant parameters explicitly. 

As noted in Sec. IIIBl the momentum space approach may become invalid if the coupling 
becomes so strong that modes far from the Fermi surface start to play an important role. 
This will be easily noticed, since it leads to A-dependence of the low-energy observables. 

We expect the effective theory (fTOj) to have a weaker sign problem than the full theory 
([2]) because it is closer to having a particle-hole symmetry: we have discarded the degrees 
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Pf 
m(A) 



P 



M- 



U 



Pf 
m{A) 



P 
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of freedom far from the Fermi surface, which, in the full theory have very different phase 
space and dispersion relations above and below the Fermi surface. As we will see below, in 
the effective theory there is a variable transformation that almost guarantees the positivity 
of the fermion determinant: it is only violated by the pu — qv piece in the interaction term. 
In Sec. mil we will discuss how this piece might affect the fermion determinant; it is formally 
subleading in p/pr, but for now we keep it in the action because we do not know a priori 
the typical scale of variation of the relevant A-modes, so it may not be rigorously negligible. 

We expect low-energy degrees of freedom to be well described by the leading-order fermion 
matrices ffTOl) . This theory still has a sign problem: det Mferm = det Mi| det Mi| is in general 
complex. To see why, consider the complex conjugate of det(Mi|): 

det (Mil)* = det(M*^) 

det((dt + ^ p)Su,vSp,g - iA*{t,pF{u - v) + pu - qv)^ 

det[{dt + ^ p)Su^ij6p^g - iA{t, -Pf{u - v) - pu + qv) 



(11) 



det ({dt-^ p')5ufi5p',q' - iA{t, -Pf{u - v) + p'u - q'v] 



where we have used the fact that the auxiliary field A is real, so A*{t,p) = A{t,—p), and 
that the relabelling p' = —p, q' = —q does not change the determinant. Compare this with 
the determinant for the spin-down particles 

det(MiJ = det (^{dt - ^ p)Su^i,6p^g - iA{t, -ppiu - v) - pu + qv)^ . (12) 

For a general A-field configuration, det(Mi|)* differs from det(Mij), so the total weight 
det(Mi|) det (Ml j) will not in general be positive semi definite. However we can see that 
det(Mi-|-)* is some sense "almost" equal to det(Mi|). Only the last term in the argument 
of the A field differs between ( !T2!) and the last line of ( fTTl) . If the offending terms could be 
ignored or suitably modified (for example if the p/pr expansion converged well enough, or 
if scattering between patches proved to be negligible — see below) then the fermion matrix 
would be positive-semidefinite. 



III. MODIFIED EFFECTIVE ACTION WITHOUT A SIGN PROBLEM 

To test the idea that the effective theory flTUl) has a tractably weak sign problem, one must 
modify it to obtain an exactly positive-semidefinite weighting, and calculate observables 
with a reweighting factor that gives results appropriate to the action flTUl) . If the sign 
problem is sufficiently weak, the reweighting factor will not be too noisy, and useful results 
can be obtained. One approach would be to generate ensembles weighted by | det Mi | or 
I Re det Mi|, and then reweight back to the proper det Mi weighting. This has not yet been 
attempted, perhaps because in the position-space formulation ]1, [6| the calculation of det Mi 
requires the construction of a non-local operator which is used in the coupling of the fermion 
modes to the auxiliary field. We suggest that formulating the theory in momentum space, 
as we have done in this paper, may make it easier to perform numerical investigations of the 
high density effective theory. For example, it is easy to do small-scale studies using a small 
number of angular patches on the Fermi surface. 

In the remainder of this section we will describe an alternative to using | det Mi | or 
I Re det Ml I, namely modifying the interaction term in ffTOl) so as to make the fermion de- 
terminant positive-semidefinite. 
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A. Action with no scattering between patches 



The simplest and most radical approach would be to make the different directions on the 
Fermi surface independent of each other, by coupling fermion modes only to auxiliary field 
modes that scatter them within the same angular patch on the Fermi surface. 

To see that such a theory is positive semidefinite, consider a theory in one spatial di- 
mension [d = 1) where the angular variable u can take on only two values, a and —a, and 
we discard high-momentum A-field modes that could scatter fermions from d to —d. This 
situation was studied in Sec. II of Ref. ^ (see also [3]). The fermion action for a single 
spin state is 



^flL= / dpdqip\d,p)(6{p-q)dt + iA{{p-q)d) + —p6{p-q)jilj{d,q) 

dpdqijj\-d,p) (6{p - q)dt + iA{{p - q){-d)) + —p6{p - g) )■?/'(-«, q) . 
\ m / 



In each line, the first two terms are anti-Hermitian and the last one is Hermitian. Now, in 
the —d patch (second line of (fT3|) ) we change variables to p' = —p and q' = —q and then 
rename these back to p and q, giving 

'^fJrm= / dpdqip\d,p)(5{p - q)dt + iA{{p - q)d) + p6{p - q))ip{d,q) 

f f \ ^''^ 

+ / dpdqil}\-d,-p){5{p-q)dt + iA{{p-q)d)-p5{p-q)\il){-d,-q) . 

Note that the fermion modes living at —a are different from the modes at a, but they couple 
to the same auxiliary field modes. Equation 0141) leads to a fermion matrix of the form 

^^■^ -[ Ma- Mh J ^ ^ 

where Mh is a Hermitian operator [p in this case) and Ma is an anti-Hermitian operator 
[dt + iA in this case). We assume that both are of even dimension. It then follows that 
M^'^^ has positive-semi definite determinant, since the eigenvalues of Ma — Mh are (up to a 
sign) the complex conjugates of the eigenvalues of Ma + Mh. The spin-down matrix M^^^ 
is independently positive-semidefinite, via a similar argument. This agrees with the result 
obtained in Ref. |6|. 

It was important in obtaining this result that the off-diagonal blocks in ( fT5|) were zero: 
in a d = 1 theory this restriction corresponds to putting a cutoff Aa on the auxiliary field 
modes, where A <^ A^ -C Pf, which imposes the condition that the A-field can only scatter 
fermions within angular patches. In higher dimensions, however, there is no separation 
of scales between momenta that scatter between patches and momenta that scatter within 
patches: an A-field of arbitrarily small momentum can scatter a fermion from one patch to a 
neighboring one, and the fermion matrices are no longer positive semidefinite (see Sec. IA 3p . 



The argument made in Ref. [6| is that, in QCD, even in more than one dimension, scattering 
between neighboring patches is subleading in fiR (where R is the in-medium screening 
distance for gluons) so scattering between patches can be ignored, and high-density QCD can 
be treated as a collection of d = 1 theories. However, in a first principles non-perturbative 
approach one cannot rely on such essentially perturbative results. The screening distance is 
a priori unknown and should be determined in a simulation. Therefore, in this article we 
will not assume that scattering between patches can be neglected. 
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B. Action that is symmetrized in the residual momentum 

We propose an alternative way of making the fermion determinant positive semidefinite: 
modify the couphng to the background field by symmetrizing it in the residual momentum 
p. This means using a new fermion determinant Mis which is constructed by replacing the 
background field A in fllOl) with 

As{p — q) = - (^A{pf{u — v) + pu — qv) + A{pf{u — v) — pu + g£>) j . (16) 

The fermion determinant with this symmetrized coupling to the auxiliary field is guar- 
anteed to be non- negative for any A(p). This follows from the argument of section 
III CI by coupling to Ag we make the interaction term symmetric in the residual mo- 
menta p and g, so now det Mi^j^ = (detMisf)*, and the symmetrized full determinant 
det Mferm,s = det Mis^ det M^i is real and positive for any A-field configuration. 

The modified ensemble generated using As has larger violation of momentum conservation 
(see Appendix which will be corrected by reweighting back to the original fermion matrix 
(see Sec. |V]). Since the modification is (formally) subleading in p/pp and hence A/pp, the 
reweighting factor may be sufficiently well-behaved that the modified theory is still suitable 
for calculations of low-energy observables in the original theory. That is one of the main 
conjectures of this paper, and we hope it can be tested by explicit Monte-Carlo calculations. 

A more radical modification of the fermion action would be to retain only the leading 
piece in the expansion, A{pp{u — v)), throwing out all dependence on the residual momenta. 
This would correspond to having all fermion modes in a given patch u couple to all modes 
in another patch v via the same auxiliary field mode. This might serve as a good starting 
point in numerical investigations because of its relative simplicity. 

It needs to be emphasized that, unfortunately, the proposed approximation may not be 
systematically improved. The approximation hinges on the observation that the theory with 
couplings symmetrized in p (which corresponds to discarding terms in its Taylor expansion 
that are odd powers of p) is positive semidefinite. Inclusion of those odd powers and other 
higher-order terms in p/pp will inevitably re-introduce the complexity existing in the un- 
modified theory ffTOl) . even more so in ([7]), ([8]). However, our conjecture is that simulations 
based on the modified theory described above will capture low-energy properties of the full 
theory with a controllable error of order of T/fi. 

IV. LATTICE REALIZATION OF THE THEORY 

To numerically evaluate expectation values of observables by Monte-Carlo methods, it 
is necessary to formulate the theory on a lattice. We will discretize the Euclidean time 
direction in the standard way, with a lattice of Nt points covering the range [0, 1/T]. The 
spatial directions will be discretized as a lattice in momentum space. As described above we 
use polar co-ordinates: we discretize the angular direction u and the radial direction p, and 
keep only modes near the Fermi surface by imposing an ultraviolet (UV) cutoff A on the 
residual momentum p. We choose a single infrared (IR) cutoff 6p, which is the momentum 
lattice spacing in the radial direction, and we choose this to also be the momentum lattice 
spacing in the angular direction on the inside edge of the momentum shell (at p = —A). 
An example of such a momentum lattice for fermions in two dimension is shown in Fig. [H 
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FIG. 1: Discretization of the space of fermionic momentum modes near a two-dimensional Fermi 
surface. This example has momentum lattice spacing (infrared cutoff) 6p = pp/8, and ultraviolet 
cutoff A = pf/4:. 

In order to simulate properties of infinite volume continuum theory the lattice parameters 
need to satisfy 

A > Ptyp > Sp , (17) 

where ptyp is the characteristic momentum scale of the physics of interest, for example the 
thermal momentum mT/pp, or a momentum carried by a correlation function. 

The modes of the auxiliary field A are most conveniently chosen to lie on a cubic lattice 
of lattice spacing 6p in all spatial directions. In order to couple them to fermion modes 
that live on a polar momentum lattice, while maintaining the correct large-volume {Sp — 0) 
limit, we must allow a small violation of conservation of momentum, which goes to zero 
a.s 6p ^ 0. This arises because, in coupling the auxiliary field to the fermions, we have 
to treat momenta that differ from each other by an amount less than half the momentum 
resolution 6p as the same momentum. Then a pair of Fermion states labeled by u, p and 
V, q corresponding to a momentum difference Ap = Pf{u — v)+pu — qv, will be coupled to 
all cubic lattice modes A(k) such that | Ap — k| < 6p/2. 

Momentum violation at the scale of an infrared cutoff is not unusual: it arises when 
one performs lattice calculations in a box with hard-wall or reflecting boundary conditions, 
in which case momentum conservation is violated at distance scales on the order of the 
size of the box, but this artifact disappears in the large- volume limit. For a more detailed 
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discussion, see appendix |X1 From (fTOj) . our discretized fermion matrices are 



Mjf = (^dt + ^ 5u,i,5p^q + I ^ A{t, k)W{pF{u -v)+pu- qv - k) 

f X ' (18) 

Mjf = f 5t - ^ P j - ^ 5^ A{t, -k)W{-VF{u -v)-pu + qv-k) . 

k 

In these expressions, W{q) (see (lAip ) is a function that imposes approximate momentum 
conservation. It is is peaked around |g| = 0, with a width ^ 6p. The role of the W 
function is to define a mapping between the A-field momenta k, which lie on a cubic lattice 
of spacing Sp, and the polar lattice momentum differences Ap = Pf{u — v) + pit — qv, such 
that \k- Ap\< 6p/2. 

As in the momentum continuum (Sec. IIIip . the lattice fermion operators flTS]) have non- 
positive determinants, but one can modify the coupling to the background field so that 
a given pair of fermion fields ■ip^{u,p)ip{v,q) couples to the symmetrized combination of 
auxiliary fields ^{A{pf{u — v) + pu — qv) + A{pf{u — v) — pii + qv)}. 

The result is a pattern of coupling between the A-field and the fermions that has the 
fermion matrices Mis, where 

^uf = + ^ P)Su,v^,<i +'^Y1 ^) ( ^iPFi^ - +pu-qv-k) 

k 

+ W{pf{u — v) — pu + qv — k) 
Mlfi = {dt- — p)Su,vSp,, -'-Y, A{t, k) ( W{-pf{u -v)-pu + qv-k) 



(19) 



+ W{—Pf{u — v) + pu — qv — k) 

It is then straightforward to verify, following the same procedure as in Eqns. (ITT]) and (1121) . 
that det(Misi)* = det(Mi5j^), so the fermion matrices (I19p have no sign problem, and their 
determinants can be used as Monte-Carlo weights to generate an ensemble of configurations. 
As noted in Appendix \^ in the modified theory an auxiliary field mode with momentum k 
couples not only to pairs of fermion modes whose momentum difference is within 6p of k, but 
also to pairs of fermion modes with momentum differences that deviate from k by as much 
as A. This additional momentum non-conservation should be corrected by the reweighting 
back to the original fermion matrix (see Sec. |V]), and it remains to be seen whether it will 
lead to a significant sign problem in computing the reweighting factor. 



V. CONCLUSIONS 

We have formulated a theory of non-relativistic fermions with a spin-independent repul- 
sive interaction in polar coordinates in momentum space. This makes it very straightforward 
for us to write down the high-density effective theory for the fermionic modes with the lowest 
free energy (those near the perturbative Fermi surface) which are expected to be relevant 
for low-energy properties of the system. We have shown that in general this effective theory 
still has a sign problem. We argued that various modifications of the action would yield 
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a positive semidefinite weighting factor Wpos which is suitable for Monte-Carlo sampling. 
These include completely decoupling different angular patches to give a set of decoupled 
l+l-dimensional theories (Sec. IIII Al and [6|); symmetrizing the interaction in the residual 
momentum p, i.e. the momentum distance from the Fermi surface (Sec. IIII B|) : and taking the 
modulus of the fermion determinant. The last two modifications are (formally) subleading 
in p/pf. 

To obtain physical results one would have to re-weight back to the original theory, by cal- 
culating the reweighting factor TZ = (det Mi/iypos)pos, where Mi is the fermion determinant 
for the unmodified effective theory of the low-energy degrees of freedom ffTOj) . and (■ ■ ■)pos is 
an average in the ensemble weighted by Wpos- If TZ turns out to have sufficiently small fluctu- 
ations on lattices whose volume is big enough to be physically relevant, then the reweighting 
procedure can be used to calculate observables in the original theory. A priori, we cannot 
say which of the modifications proposed here will yield a positive ensemble with the closest 
overlap with the physical ensemble (i.e. the smallest fluctuations in TV), or whether any of 
them will be close enough to be useful. This is a typical issue with reweighting schemes. One 
may get some information by comparing simulations that use different modifications. Our 
approach reduces the original theory to a Yukawa-type theory of the degrees of freedom near 
the Fermi surface, and it is encouraging that such theories have been found to have a fairly 
gentle sign problem, so that simulating using the modulus of the fermion determinant did 



not introduce unphysical phases [15|]. Rough preliminary investigations of a two-dimensional 
model (see Appendix [B]) seem to indicate that in the momentum space formulation the sign 
problem is indeed very moderate. The momentum space approach may become invalid if 
the coupling becomes so strong that modes far from the Fermi surface start to play an 
important role. As we noted in Sec. Ill Bl it will be straightforward to check whether this 
is a problem, since such a breakdown will be signaled by A-dependence of the low-energy 
physical observables. 

Working in momentum space gives the theory a perfect kinetic term. We keep only the 
leading order in p/pr but there are no discretization errors coming from approximating 
derivatives by finite spatial differences. The polar momentum formulation also makes it 
straightforward to develop simplified models for testing the sign problem in the effective 
theory, such as ones with an auxiliary field that takes on discrete values (see Appendix [B]) 
or with a small number of patches on the Fermi surface (see Sec. lA 2\i . However, it imposes 
some costs. Firstly, as discussed in Appendix \^ we have to introduce an infrared cutoff Sp, 
and in the interaction between the auxiliary field and the fermions we treat fermion modes 
whose momenta that are within ^Sp of each other as if they had the same momentum. This 
means the interactions introduce momentum violation of order 6p. Secondly, our formulation 
is more expensive for lattice computations because the interaction term, which was local 
in position space, becomes non-local in momentum space. This means that Mi and any 
modified version of it are non-sparse matrices. (Note, however, that the position-space 
formulation j^, [6| also requires the calculation of a non-local operator). Computation of 
the determinant of a generic non-sparse N x N matrix takes of order A^'^ operations [l7| . 
whereas the algorithms for dynamical fermions with local actions typically take of order N 
operations at fixed UV cutoff [l8|. This extra cost will easily be worthwhile, however, if the 
resultant ensemble has a less severe sign problem. The number of configurations required for 
reweighting from a positive action to the original action rises extremely quickly, as exp(A^A/) 
(where A/ is the free energy density difference between the actions |3|), so if Wpos has a 
smaller value of A/ than existing approaches then this could easily compensate for the extra 
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costs arising from the non-sparse fermion matrix. We hope that lattice calculations using 
the high-density effective theory will test these ideas in the future. 
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APPENDIX A: DISCRETIZATION IN POLAR MOMENTUM SPACE 



In this appendix we discuss some of the technical issues that arise when we discretize the 
fermion degrees of freedom in polar coordinates in momentum space. The main question 
is how the fermions should be coupled to the auxiliary field, which is also discretized in 
momentum space. 

1. Entire Fermi surface 



FIG. 2: Momentum-conserving coupling of auxiliary field A to modes near the Fermi surface, in 
two spatial dimensions. Panel (a): continuous momentum space, where an j4-mode with given 
momentum will couple to any pair of fermions with the appropriate momentum difference (thick 
lines with embedded arrow heads, red online; just a few examples are shown). Panel (b): dis- 
cretized momentum space for fermions using polar co-ordinates. If we require exact momentum 
conservation, then a given (non-radial) yl-mode can couple to only one pair of fermion modes and 
their antipodal images. 

If we discretize the whole Fermi surface then we cannot maintain exact momentum con- 
servation in the interactions that are mediated by the auxiliary A field. We illustrate the 
issue for the two-dimensional case in Fig. [2j In the infinite volume limit (continuum in mo- 
mentum space) the A-field mode with a given momentum k can couple to any two fermion 
modes whose momenta differ by k: the Lagrangian will contain many terms of the form 
ip^{p)A{k)ip{p — k) for different values of p. These pairs of modes densely cover the available 
momentum space, as shown in Fig. [2^ by the many ways one can draw lines of momentum k 
(thick lines with embedded arrow heads) connecting pairs of fermion modes near the Fermi 
surface. However, once one discretizes the radial and angular components of the momenta, 
the momentum difference between any two fermion modes becomes unique (up to inver- 
sion through the origin, and apart from purely radial momentum differences). So in the 




(a) 



(b) 
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FIG. 3: Large-volume {5p — > 0) limit showing which pairs of fermion modes couple to an ^-mode 
with momentum (kx, ky) = (3, 2), if we tolerate momentum non-conservation ca = ^^P- We work 
in units where pp = 10 and A = 4. The four graphs correspond to 5p = 2.67, 5p = 1.6, 6p = 1.14, 
and 6p = 0.889. 

discretized theory (Fig. ^jp) if we require exact conservation of momentum then a typical 
y4-mode either cannot couple to any fermion modes, or can only couple to one pair of fermion 
modes, and their antipodal images. This remains true no matter how small Sp becomes. 
We conclude that an exactly-momentum-conserving discretization will not give the correct 
theory in the 6p —>■ limit. For example, if one made the coupling constant g small, most 
fermion-fermion scattering amplitudes would be proportional to g in the full theory (from 
single A exchange) but would have no 0{g) contribution in the discretized theory, where 
single A exchange would be forbidden by momentum conservation. 

We suggest that this problem can be solved by allowing a small tolerance for mo- 
mentum mismatch: we will allow a given A mode of momentum k to couple to any 
pair of fermionic modes whose momenta differ by an amount Ap that is within of k, 
i.e. I Ap — k\ < eA- For simplicity, we will choose the A-modes to live on a cubic lattice with 
the same momentum lattice spacing 6p as the fermionic lattice, and we will set the tolerance 
eA = ^Sp, so that typically a pair of fermionic modes will couple to at most one A-mode. 
The yl-mode momentum lattice will have a UV cutoff of 2(pi7' + A), so that fermion modes 
can be scattered between distant parts of the Fermi surface. We expect that in the limit 
5p — > 0, infrared artifacts (arising from the momentum lattice spacing and the tolerance 
for momentum non-conservation) will disappear. In Fig. [3] we show how this scheme solves 
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(a) 



(b) 



FIG. 4: Figure showing which pairs of fermion modes couple to a the A-field mode with momentum 
{kx,ky) = (3,2), for pp = 10, UV cutoff A = 4, IR cutoff 6p = 2. Panel (a): unmodified theory 
(see (|18p ): momentum violation is of order 6p. Panel (b): theory with interaction symmetrized in 
p (see (|16p . (|19p ) to ensure a positive semidefinite fermion determinant; momentum violation is of 
order A. 



the problem illustrated in Fig. [2l In Fig. [3|, as 5p — 0, a given A-mode, in this case the 
one with momentum {kx,ky) = (3,2), couples to more and more pairs of fermion modes, 
approximating the dense set of Fig. [2^. 

In this approach, then, the discretized coupling of the A field to the spin-up fermions is 
given by 

M\t,,u,, = Yl ^)W{pf{u - v) +pu-qv-k) (M) 

k 

where W imposes momentum conservation with a tolerance of ^6p, 

{1, for Isl < 7^6p 
U, tor |s| > -^op 

In section HVl we discussed how the coupling of the fermions to the auxiliary field could be 
modified (symmetrized in the residual momentum p) to yield a positive semidefinite fermion 
matrix. In Fig. H] we illustrate this modification for one of the auxiliary field modes. We see 
that in the modified theory (Fig. Hb) the A-mode with momentum [kx, ky) = (3, 2) couples 
not only to pairs of fermion modes whose momentum difference is within about Sp of k, but 
also to pairs of fermion modes whose momentum difference deviates from k by up to ±A. 
This additional momentum non-conservation should be corrected by the reweighting back 
to the original fermion matrix, and it remains to be seen whether it will lead to a significant 
sign problem in computing the reweighting factor. 
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(a) (b) 

FIG. 5: Momentum-space discretization of modes in two antipodal small patches on the Fermi 
surface. Panel (a) shows the modes in polar co-ordinates; panel (b) shows the same set of modes 
realized as a flattened Fermi surface, in which m is a periodic Cartesian coordinate orthogonal to p. 
The double-headed arrows (red online) show some of the couplings of the A-modes with momentum 
{k^,ky) = {3,2)6p. 



2. Small antipodal patches 

A less complicated coupling of the A-modes to the fermions is possible if we restrict 
ourselves to a pair of small antipodal patches on the Fermi surface. In Fig. [5t^a) we show 
how this would work in two spatial dimensions, where p gives the radial momentum relative 
to the perturbative Fermi surface, and u is the polar angle of the momentum. 

The approximation of Eq. ( fTOl) neglects the difference between the amount of phase space 
above and below the Fermi surface, and can therefore be envisaged as a flattening of the 
Fermi surface to a line, as shown in Fig. MJo), where u becomes a Cartesian momentum 
co-ordinate orthogonal to the radial momentum p. In each patch we treat u as a periodic 
variable. For simplicity we assume that the lattice spacing has the same value 6p in both 
angular and radial directions, but one could use an anisotropic lattice if necessary. The 
generalization to three dimensions is straightforward: simply add another angular variable. 

In Eq. ( 1101) . each mode of the A-field couples equally to all pairs of fermion modes whose 
momenta differ by exactly the momentum of the A field. We implement this in the lattice 
theory after the flattening of the Fermi surface, so the A field lives on a similar momentum 
lattice, with the same lattice spacings, generated by taking all possible momentum differences 
on the lattice of Fig. MJo)- Following Eq. (fTUi) . we allow each mode of the A-field to couple 
to all pairs of fermion modes with the right momentum difference. In Fig. El^b) we show 
some of the pairs that one specific A-mode would couple to (double-headed arrows). When 
one takes in to account the periodic boundary condition in the u direction, these all have 
the same momentum difference. 

The small patch scheme is a radical truncation of the full theory, because it only keeps 
a small part of the Fermi surface. One could regard it as a crude way of implementing a 
forward-scattering-dominated interaction. The ground state of a single-patch theory is more 



likely to be a particle-hole condensate [19|, |20| than a condensate of Cooper pairs. But it 
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could also be used as a tractable toy model for testing the ideas that we have outlined in this 
article. It is simpler than the "whole Fermi surface" scheme of the previous section because 
there is a smaller number of degrees of freedom, and momentum is exactly conserved. 



3. Positivity tests 




FIG. 6: Simplified discretization including only fermion modes with momenta i{pF ^ p)a and 
i{pF ip)b. These couple to six auxiliary field modes with momenta ki . . . kg. We have shown ^3 
and ^4 in the antipodal patches to minimize clutter. 



We can test the positivity of the low energy effective theory in a very simple example, 
such as that depicted in Fig. [61 We focus on one spin component of the field, and reduce 
the fermionic field to two patches, a and b, and their antipodal partners. In each patch we 
have two residual momenta ±p, spaced equally above and below the Fermi surface. Even if 
we decouple high-momentum A-field modes, so that there is no coupling between the (a, b) 
patches and the (—a, —b) patches, the fermion determinant is still not positive definite. To 
see this, note that the fermion matrix is block diagonal 



M+ 
M_ 



(A3) 



where M+ is the fermion matrix for the (a, b) patches and M_ = (M+)^ (see below) is the 
fermion matrix for the {—a, —b) patches. The full fermion determinant is therefore 



det(M) = det(M(^)) det(M(^)) = det{M^^^y , 
In the basis {il){a, +p), ?/'(a, -p), ^^(6, +p), ■il){b, -p)). 



det(M(^)) = det(M+)^ 



(A4) 



P + ik 



( p %A{kx) 
%A[—k]) —p 
iA{-k2) iA{-k3) 



%A{k2) 
iA{k^ 
p 



\iA{-k3) %A{-k^) tA{-ke 



zA{h) \ 
iAih) 
iA{ke) 
-p ) 



(A5) 



where, as shown in Fig. O /ci = —2pa, k2 = {pr + p){b — a), etc. All the ki are much less 
than Pp. The residual momentum term, P, is a diagonal Hermitian matrix with entries ±p, 
and A is an off-diagonal Hermitian matrix because A{—k) = A{k)*. We can see that M_ = 
P -|- zA* = because A field momenta are reversed when we perform the transformation 
to the antipodal momenta. It is easy to verify that for a generic A-field configuration, 
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det(M+) is a generic complex number, so the full fermion determinant (lA4p of the effective 
theory is not real. If we turn off the coupling between neighboring patches by setting 
A{k2) = ^(fca) = ^(^4) = ^(fcs) = then the determinant is always real, in agreement 
with the argument of Sec. IIII A[ 



APPENDIX B: TESTING USING A TWO-DIMENSIONAL EXAMPLE 



In this appendix we describe a simple theory that can be used to test the ideas set forth in 
this paper. We start with fermions in two spatial dimensions, feeling a point-like repulsion 



Hint = g V'tW^TWV'i(x)^/'i(x) 



(Bl) 



where g > 0. We employ Hirsch's variant of the Hubbard- Stratonovich transformation, 
where the auxiliary field a is discrete, taking on values ±1. From the start we assume the 
Euclidean time direction is periodic with period T, and discretized with time spacing r. 
Then we find [131 



exp[-r^nt(f)n|(f)] = ^ ^ exp[26a(f)(n^(f) - - y + (B2) 

a{x)=±l 

where na{x) = ■ipl^{x)ipa{x), and a = {t, i}, and 



b = — , tanh^(6) = tanh ( — 
T V 4 



(B3) 



The fermion matrices in the effective theory in momentum space are then given by 



M 



Ft 



M 



Ot^ P 

m . 

Ot p 

m . 



5p,q5u,v + 2ba(t,pF{u - v) + pu - qv), 
5p,q5u,v + 2ba{t, -Pf{u - v) - pu + qv), 



(B4) 



where Pf = y^2mjl and ft = n — g/2 as follows from flB2l) . Note that we use a Nambu- 
Gor'kov-like basis where the a{t, p — q) term is positive in both the spin-up and spin-down 
fermion matrices. The auxiliary field does not have any action of its own {Saux = 0). 

We now discretize the system for lattice calculations. The radial residual momentum has 
an infra-red cutoff 6p and ultra-violet cutoff A, and the angular variable is discretized in 
steps (see Sec. lIVp so the momentum lattice action is 



PfL 



k,n=—Nr, u,v=0 



Pk 



PFPn 



m 



~PfL 

2 b——a[pF[u -v)+ pkU - pnv) 



(B5) 



where Np = A/6p, pn = n6p, v is a unit vector at angle 2'Kv/Ny to the x-axis, we have 
defined the IR cutoff length scale L = 27T/6p, and Vj is a discretized version of the time 
derivative. The fermion matrix dimensionality is Nt x Ny x {2Np + 1). 

The auxiliary fields, a{t, p), are defined on a square lattice in Cartesian momentum space, 
with UV cutoffs ±Ao- on the x and y components of momentum, where A^. = -|- A, and a 
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momentum lattice spacing (IR cutoff) (5po-, and we define L^j = 27r /Spa- In practice we work 
with a position space square lattice on which cr(t,x) = ±1 and we Fourier transform it to 
obtain a(t, p). We set 6pa = Sp. The number of sites along each side of the a-lattice is an 
odd integer, 2A„ /Sp^j + 1. 

We have performed very preliminary calculations for pp = 2.58, m = 1.0, /i = 3.33, 
T = 1.0, with coupling g = 1.0. For the fermions we used a momentum space lattice with 
6p = 0.194, A = 0.775, = 64. This corresponds to Np = 4, A/pp = 0.3, T/fi = 0.3. 
In the Euclidean time direction we set r = 0.111 which corresponds to Nt = 9. The 
dimension of the fermion matrix is 5184. For the auxiliary field we use A^r = Pf + A = 3.35 
Spa = Sp = 0.194, so the a field lives on a 35 x 35 cubic lattice. 

To get an idea of the severity of the sign problem in this formulation of the theory, we 
performed a Monte-Carlo calculation of the reweighting factor TZ (see Sec.|V]), We formulated 
the theory in Mathematica using its Det function to numerically evaluate the determinant of 
the 5184 X 5184 fermion matrix. We used the Metropolis algorithm, with an accept/reject 
step for each update of the auxiliary field a using the non-negative weighting factor Wpos = 
I det MiTi-MiTjl. We note that even with this numerically unsophisticated approach, we are 
able to study a lattice that covers an entire 2D Fermi surface, with reasonable resolution (A 
four times bigger than 6p (ITTl) ). 

The observable that we measured is the real part of the reweighting factor, ReT^ = 
cos(argdet Mir|Mi?j^(o")). In our Monte-Carlo runs we started with all the cr(t, x) being 
-|-1, and found that when we had done enough updates (about 10"^) the a were equally 
distributed between -|-1 and —1, and the average value of Re7?.(cr) — 1 was of order 10~^, 
with occasional fluctuations no larger than about 10~^. These initial results have not been 
checked for proper thermalization, autocorrelation, and finite volume effects, but they are an 
encouraging indication that the theory of the degrees of freedom near the Fermi surface has 
a moderate sign problem, with a reweighting factor that is close to 1, rather than fluctuating 
around zero. 
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